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Using the time-dependent density-matrix renormalization group (tDMRG), we study the time evolution of 
electron wave packets in one-dimensional (ID) metal-superconductor heterostructures. The results show An¬ 
dreev reflection at the interface, as expected. By combining these results with the well-known single-spin- 
species electron-hole transformation in the Hubbard model, we predict an analogous spin Andreev reflection in 
metal-Mott insulator heterostructures. This effect is numerically confirmed using ID tDMRG, but it is expected 
to be present also in higher dimensions, as well as in more general Hamiltonians. We present an intuitive picture 
of the spin reflection, analogous to that of Andreev reflection at metal-superconductors interfaces. This allows 
us to discuss a novel antiferromagnetic proximity effect. Possible experimental realizations are discussed. 

PACS numbers: 71.10.Fd, 71.10.Pm, 74.20.-z, 74.45.+C 


Introduction .—Correlated electrons at the interface of two 
materials can exhibit a wide range of remarkable phenom¬ 
ena. Among the most interesting effects is the Andreev reflec¬ 
tion (AR) at the normal metal-superconductor (N-SC) inter¬ 
face [1]. As an electron is transmitted from the normal metal 
(N) into the superconductor (SC), it attracts a second electron 
from the metal to form a Cooper pair. The second paired elec¬ 
tron leaves a hole that is reflected off the interface into the 
metal. In AR a charge current in the metal becomes a super¬ 
current as it enters the SC. 

Several variants of AR have been widely discussed in recent 
years. Specular AR has been predicted at graphene N-SC in¬ 
terfaces, where the reflection angle is inverted [2]. Proximity 
effect and AR in SC-ferromagnet heterostructures have been 
studied both theoretically [3, 4] and experimentally [5-7]. 
Spin-dependent Q-reflection was predicted in normal metal- 
itinerant antiferromagnet (N-AFM) interfaces [8], where a tt 
phase shift is seen between the reflected spin-up and spin- 
down electrons. In N-SC-N systems, nonlocal or crossed AR 
(CAR) involves the transmission of an electron from a metal 
to the SC and the creation of a hole in the other metal. In trans¬ 
port experiments, CAR competes with elastic cotunneling as 
the dominant mechanism [9]. Due to its nonlocal nature, CAR 
can generate nonlocal entanglement and correlations in car¬ 
bon nanotubes [10] and topological superconductors [11]. 

In one-dimensional (ID) systems, dominant superconduct¬ 
ing fluctuations are characterized by the Luttinger liquid (LL) 
parameter Kp > 1, or AT > 1 in the case of spinless 
fermions [12]. AR was predicted for ID spinless fermions 
with nearest-neighbor (NN) attraction U < 0 coupled to a 
noninteracting wire {K = 1). The reflection coefficient is 
given by 7 = {l — K)/(lf-K). ¥or K > 1, 7 < 0, which cor¬ 
responds to AR [13]. This formula can be generalized to the 
spin and charge sectors of electrons in ID [14]. The spinless 
fermions results were later confirmed using time-dependent 
density matrix renormalization group (tDMRG) calculations 


that cast the AR in the context of ID cold atoms [15]. To our 
knowledge, a real-time study of AR in the case of spinful elec¬ 
trons has not been presented before. We are also unaware of 
any prediction of its spin analogue described in this work. 

In this Letter, we propose a novel spin Andreev-like re¬ 
flection (SAR) in normal metal-Mott insulator (N-MI) het¬ 
erostructures. An electron with a given spin projection under¬ 
goes a spin flip upon reflection, inducing a spin-1 excitation 
in the MI. We present an intuitive picture of the process, anal¬ 
ogous to that of AR in N-SC. SAR is verified in a Hubbard 
chain using tDMRG. However, the effect is not restricted to 
ID, and can be observed in higher dimensions as well. We 
discuss the differences of our work with previous efforts, and 
present possible experimental realizations of the new effect. 

The organization of our work is the following. First, we 
study the dynamics of a wave packet with both spin and charge 
colliding with a N-SC interface. The superconductor is mod¬ 
eled using the Hubbard Hamiltonian with onsite attraction 
U <0. Away from half-filling, we observe the partial AR 
of the charge component. Due to the spin gap, the spin of the 
wave packet undergoes normal reflection. Then, by using a 
well-known single-spin-species electron-hole transformation, 
these results are translated into the novel SAR in N-MI sys¬ 
tems. Finally, this translation is confirmed by calculating the 
time evolution of the wave packet colliding with the N-MI in¬ 
terface and clearly observing the spin Andreev reflection. We 
explore an AFM analogue of the SC proximity effect in the 
N-MI heterostructure. To close, we propose experimental re¬ 
alizations to test our predictions. 

Model .—Using tDMRG [16-18], we study a Hubbard 
chain of Lj sites connected to a non-interacting lead of 
sites. The Hamiltonian can be written as 77 = Hj + Ht, 
where Hi represents the interactions in the Hubbard chain 
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FIG. 1. Snapshots of the charge and spin propagation for U = —4:th 
in the doped case. The spin behavior is the same as in the half-filled 
case since the spin sector is still gapped, but the charge behavior is 
different. As it enters the interacting region, the charge of the wave 
packet increases, and consequently a small wave packet of opposite 
charge is formed in the lead, which is a clear signature of partial AR. 
This can be seen between t — 12 and t = 20. The reflected charge 
and spin move at the Fermi velocity of the lead. 

and Ht is the kinetic energy of the entire system, 

HT = -th (cU+i.+H.c.). (2) 

cr,2=l 

The hopping integral th is taken as the energy unit. The chem¬ 
ical potential /i of the Hubbard chain relative to the lead con¬ 
trols its charge density. The rest of the notation is standard. 
Note that H has SU(2) symmetry in the spin sector and, for 
/i = 0, the same symmetry in the charge sector. For details of 
the implementation of the time evolution of the wave packets 
see Refs. 19 and 20. 

Results. —For U < 0, we study the half-filled case with 
/i = 0, and = 30, and a doped case with 

= 24 and /i = —0.2. At half-filling, the charge 
is perfectly transmitted and the spin is totally reflected due to 
the absence and presence of energy gaps, respectively. There¬ 
fore, we observe no evidence of AR. This can be understood 
by noticing that, for U <0, Kp = 1 (7 = 0): the free-fermion 
limit [20]. 

The results for the doped case and U = — are shown in 
Fig. 1. The behavior of the spin component is the same as in 
the half-filled case, since the spin sector is still gapped. The 
charge behavior is more interesting. After entering the SC re¬ 
gion, the charge of the wave packet increases, while a small 
wave packet of opposite charge is formed at the lead. This is a 
clear signature of partial AR, i.e. the reflection of an opposite 
charge at the interface. This reflected charge and spin move at 
the Fermi velocity of the lead. Note that in this case Kp > 1, 
leading to a negative reflection coefficient: this means that ad¬ 
ditional charge is attracted into the superconductor, or equiv- 


FIG. 2. Wave packet propagation for the repulsive Hubbard model 
with a finite magnetization (see text for parameters used). The wave 
packet spin is transmitted, whereas the charge is refiected due to the 
charge gap. SAR is observed. As the spin enters the MI, it expels an 
opposite spin back to the lead, i.e. an effective spin flip takes place. 
The spin fiip here is partial, as in the charge AR case. This refiection 
of an opposite spin is the signature of the proposed SAR. 

alently, opposite charge is reflected. For larger values of \U\, 
such SiS U = —Sth, AR is not observed within our resolution 
and the charge undergoes normal reflection. This is consis¬ 
tent with previous studies in the regime where the SC pair¬ 
ing energy is comparable to the Fermi energy, as in ultracold 
fermionic mixtures, where it has been observed that the spec¬ 
ular reflection of particles has a robust amplitude [27]. 

Consider now the main results of our publication. A well- 
known property of the Hubbard model with onsite interaction 
is the equivalence between the charge and spin sectors of the 
U > 0 and U < 0 cases. In other words, if the Hubbard 
model with U < 0 is studied, then the U >0 properties can 
be deduced by merely exchanging charge and spin. More con¬ 
cretely, consider the electron-hole transformation on one spin 
species T: c\^ and c\^ (— 1 )*C 4 , and consequently, 

Ui^ rii^ and ^ 1 — This transformation leaves 
Ht invariant and maps Hi into H / such that 

Hi = U^K-1) (3) 

where U = —U, Sf = — n^), and B = 2/i is the 

Zeeman field. In other words, T maps the attractive Hubbard 
model into its repulsive equivalent with the charge and spin 
sectors exchanged. InHi, B breaks the spin SU(2) symmetry, 
similarly as /i breaks the charge SU(2) symmetry in Hj. 

By applying this transformation to the results shown in 
Fig. 1, or to the AR effect in general in any dimension, we 
can predict an interesting effect: an electron incident on a 
N-MI interface can undergo a spin flip upon reflection. To 
confirm this effect, the time evolution of wave packets was 
studied in the Hamiltonian H = Ht Hj. The parameters 
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FIG. 3. Charge and spin in the repulsive Hubbard model case for 
different values of U. The results are shown at t = 18. SAR is 
maximal around U = Ath- For larger U, the SAR is suppressed and 
normal spin reflection becomes dominant. Note that for small U, part 
of the charge is transmitted due to the small charge gap and the flnite 
energy of the wave packet. 


are set to those obtained from transforming H using T. That 
is, U = B = —0.4, = 24 and = 36. The SC spin 

gap translates to the charge gap of the MI. The tDMRG results 
are shown in Fig. 2. The spin behavior is identical to that of 
the charge shown in Fig. 1, and vice versa. As the spin enters 
the MI, it expels an opposite spin back to the lead. In other 
words, an effective spin flip takes place (together with a nor¬ 
mal reflection without spin flip). The reflection of an opposite 
spin is the signature of the proposed SAR. 

The effect of electron interaction on SAR is shown in Fig. 3. 
Snapshots of spin and charge propagation at f = 18 are 
shown for different values of U. In these calculations, we 
set ko = — 2a k, where kp^ is the Fermi momentum of 

the spin-up electrons, and tune B so that the magnetization in 
the system is uniform. SAR shows nonmonotonic behavior 
as a function of U. It is weak for small U/t^ < 1, where 
AFM fluctuations are unfavorable and free-fermion-like be¬ 
havior is expected. Upon further increasing U, SAR reaches 
its optimal value at intermediate coupling around U = 

For U/th ^ 1, exchange and hence long-range order are sup¬ 
pressed as J Also, normal spin reflection is observed 

in addition to SAR, and as U increases, the normal reflec¬ 
tion becomes dominant. This resulting trend is the AFM ana¬ 
logue of the BCS-BEC crossover observed in unitary Fermi 
gases [28, 29]. 

Figure 4 shows an intuitive picture of the proposed SAR 
and a comparison with the AR phenomenon. In the latter, an 
electron, with energy below the SC gap, incident on the in¬ 
terface from the metal side can be either reflected (normal 
reflection), or transmitted via a two-electron process (AR). 
The electron enters the SC in a higher energy intermediate 
state, then attracts a second electron from the metal to form a 



FIG. 4. Analogy between the AR in N-SC and SAR in N-MI. (a)-(c) 
AR mechanism: An electron, with energy below the SC gap, in¬ 
cident on the interface from the metal can be either reflected (nor¬ 
mal reflection), or transmitted via a two-electron process (AR). The 
electron enters the SC to form a higher energy intermediate state, it 
then attracts a second electron from the metal to form a Cooper pair, 
leading to the creation and reflection of a hole with opposite spin. 
(d)-(f): Analogous SAR mechanism. An electron can be reflected 
without spin exchange, or it can undergo SAR through a second or¬ 
der process. An intermediate doubly occupied state is formed, then 
an electron of opposite spin is reflected back to the metal. Note that 
the resulting net effect of AR is the transfer of two electrons into 
the SC, whereas that of the SAR is the spin increase in the MI by 
= 1. (g) Experimental setup proposed to detect crossed SAR. 
Al, A2: antiferromagnets; Nl, N2: metallic leads; d: coherence 
length; U : applied voltage; /: flowing current; V : measured voltage. 


Cooper pair, eventually leading to the creation and reflection 
of a hole with opposite spin. In the SAR case, an electron in¬ 
cident on the interface can either be reflected without spin flip, 
or it undergoes SAR through a second order process. An inter¬ 
mediate doubly occupied state is formed, and then an electron 
of opposite spin is reflected back to the metal (with a probabil¬ 
ity that depends on parameters like U). Note that the net result 
of AR is the transfer of two electrons into the SC, whereas that 
of the SAR, is the spin increase in the MI by AS^ = 1. 

A few comments are in order. (7) In Refs. 13 and 15, the 
traditional AR is predicted for ID spinless fermions with 
NN attraction V. Using the Jordan-Wigner transformation, 
this model can be mapped into a spin-^ Heisenberg model. 
The non-interacting chain is mapped into an XY spin chain, 
whereas the superconductor is mapped into a ferromagnetic 
XXZ chain. Charge AR is thus translated into a spin reflec¬ 
tion caused by the NN attraction between parallel spins. This 
provides another model in which SAR can be observed, al¬ 
though the mechanism of the reflection is different from that 
of focus in this work. (2) Reference 8 proposes an interesting 
Q-reflection at N-AFM interfaces. A quasiparticle with mo- 
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mentum k is reflected to one with k + Q where Q is the order¬ 
ing vector of the AFM. The scattering phases of spin-up and 
spin-down electrons differ by tt. However, the SAR with spin 
flip is not found mainly because of the mean-fleld treatment 
of the AFM. (3) The SAR effect proposed here can also occur 
in a three-terminal N-MI-N geometry (normal-Mott-normal), 
by analogy with the N-SC-N (normal-SC-normal) investiga¬ 
tions already mentioned [9, 30]. If the width of the SC region 
is smaller or comparable to the coherence length ^ 5 , the AR 
induced hole can be generated in the second metal leading to 
a nonlocal charge transfer. In the SAR case, a spin-up elec¬ 
tron incident over a thin Mott insulator can transmute into a 
spin-down electron on the other side. These two spins will 
be entangled, with possible applications in quantum comput¬ 
ing, as in the AR case [10, 11]. (4) Another area of potential 
value of the SAR effect is in spectroscopy where the sym¬ 
metry of the superconducting state has been studied in, e.g., 
iron-based [31] and heavy-fermion [32] superconductors via 
the canonical AR. Exotic Mott states such as d-wave insula¬ 
tors [33] could be analyzed by this procedure. In general, any 
realization currently known of the standard AR will admit a 
translation into the SAR language. 

Experimental predictions .—In recent magnetoconductance 
experiments on bilayer Aims of copper (Cu), a normal metal, 
and copper monoxide, an antiferromagnet, anomalies were 
observed when compared with Cu grown on a band insula¬ 
tor [34]. A proximity effect of antiferromagnetism inside the 
Cu layer was invoked to explain the results. This effect was 
theoretically observed before in a real-space dynamical mean- 
fleld theory computational study [35], but a simple explana¬ 
tion was not provided. In the context of the SAR effect pro¬ 
posed here, the existence of an “antiferromagnetic proximity 
effect” is natural since the canonical AR is the basis to under¬ 
stand the standard proximity effect in N-SC interfaces [36]. 
The mapping between the f/ >0 and f/ < 0 Hubbard 
model provides a simple explanation of the results reported 
in Ref. 34, and it is a concrete prediction of our effort: an¬ 
tiferromagnetic proximity effects into normal metals should 
be as ubiquitous as in the case of superconductors. Our com¬ 
putational observation that intermediate \U\ is more optimal 
than large \U\to observe the SAR lead us to believe that spin- 
density-wave (weak coupling) AFM, such as in the recently 
much investigated iron superconductors, are better than local- 
spin (strong coupling) AFM materials to test this prediction. 

An even more exotic prediction also emerges from the 
U < 0 io U >0 mapping. In superconductors the CAR, 
involving a N-SC-N interface, has attracted considerable at¬ 
tention in quantum entanglement [37]. The prediction is that 
a spin-up electron incident on the SC from the first lead may 
induce a Cooper pair by borrowing a spin-down electron from 
the second lead, as long as the superconducting width is com¬ 
parable or smaller than the coherence length. In the canonical 
CAR a hole current and an associated voltage are produced in 
the second lead. This voltage is nonlocal in the sense of be¬ 
ing in a region without a drive current [30]. Our prediction is 
that in a similar geometry, the injection of spin-up electrons 


in an AFM should produce an observable current and voltage 
of spin-down electrons in the second lead, defining a “crossed 
spin Andreev reflection” (CSAR). The four-terminal geome¬ 
try employed in Ref. 9, replacing the SC by an AFM as in 
Fig. 4 (g), provides a suitable setup to test our prediction. 

Conclusion .—In this Letter, we have predicted a novel spin 
analogue of Andreev reflection in metal-Mott insulator het¬ 
erostructures. An electron incident on the interface from the 
metal can undergo a spin flip upon reflection, thus creating a 
spin-1 excitation in the Mott insulator. This effect was verified 
in ID models using tDMRG; however, the proposed mecha¬ 
nism is valid in any dimension. Note also that the use of the 
one-spin-species electron-hole transformation in the Hubbard 
model merely provides a rapid path to the SAR notion, except 
for the case of a bipartite lattice where the transformation is 
exact. However, by mere continuity we believe that models 
that describe real materials with a density of states close to 
the particle-hole symmetric point should still display the SAR 
effect. On the other hand, if the deviation from particle-hole 
symmetry is large, for example as the doping away from half- 
filling increases substantially, then the SAR prediction must 
be revisited. Also, because of the similarity with the Andreev 
reflection, this effect can be studied in equivalent experimen¬ 
tal setups. An interesting example is using the spin reflection 
proposed here in order to generate nonlocal entanglement as 
in the crossed Andreev reflection phenomenon. 
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SI. NUMERICAL TECNIQUE 

The ground state |T^o) is calculated using static 
DMRG^’^. Then, a hole Gaussian wave packet with spin 
(7 and mean momentum is created in the lead by applying 
the operator 

^ (SI) 

3 

to l^o)^- The coefficients aj are chosen to generate a Gaus¬ 
sian wave packet centered at jo with width dx- A is a nor¬ 
malization constant. In momentum space, the wave packet is 
centered at ko with width dk = I/ctx- We set ko = kp — Adk 
and dk = 0.05, where kp is the Fermi momentum of the lead. 

The excited state 

is time evolved under H: | (t)) = I ) • We follow the 
time evolution of the charge and spin densities, 

ni{t) = niy{t) 

and 

Si{t) = i [riifit) - Unit)], 
respectively with 

so that, in the results shown, the ground state expectation val¬ 
ues are subtracted from the time-dependent measurements. 
The algorithm used in the time-dependent part of the DMRG 
is based on the Krylov space decomposition^"^. We use 350 
DMRG states in the time evolution calculations, keeping the 


relative error in the wave vector norm of the order of 10“^ 
and smaller. The calculations are done for a fixed number of 
electrons in the system, and pt is tuned so that the density is 
uniform. We study systems with Lj = Lp = 30 sites. 

S2. HALF-FILLED CASE RESULTS 

For U < 0, the half-filled case was explored for // = 0, and 
Ny = Ny = 30. Figure SI shows the results for the former. 
At t = 0, the wave packet, with charge and spin, is created 
in the non-interacting lead. The charge is transmitted without 
any refiection. The spin, on the other hand, is reflected com¬ 
pletely due to the spin gap, although a small fraction of the 
spin tunnels into the interacting region before it is reflected. 
In the case of half filling, and for U <0, the LL parame¬ 
ter is Kp = 1 which is the same as in the non-interacting 
case. Thus, the charge is transmitted without any scattering, 
and perfect transmission is expected. Similar results are also 
obtained for larger \U\. In this case, the ground state of the 
Hubbard model consists of alternating empty and doubly oc¬ 
cupied sites. As the charge enters, it moves with an effective 
hopping 



FIG. SI. Wave packet propagation for U = —Ath at half filling, rii 
and Si are shown at different times. At t = 0, the wave packet, with 
charge and spin, is created in the non-interacting lead. It reaches 
the interface around t = 8. The charge is transmitted without any 
reflection, but the spin, on the other hand, is reflected completely due 
to the spin gap. 
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